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Abstract 



We present the results of a calculation of differential cross sections and po- 
larization observables for proton-deuteron elastic scattering, for proton labo- 
ratory energies from 2.5 to 22.5 MeV. The Paris potential parametrisation of 
the nuclear force is used. As solution method for the charged-composite par- 
ticle equations the 'screening and renormalisation approach' is adopted which 
allows to correctly take into account the Coulomb repulsion between the two 
protons. Comparison is made with the precise experimental data of Sagara et 
al. [Phys. Rev. C 50, 576 (1994)] and of Sperison et al. [Nucl. Phys. A422, 
81 (1984)]. 
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I. INTRODUCTION 



Calculation of proton-deuteron (pd) scattering represents one of the most challenging 
remaining tasks in few-body nuclear physics. The interest arises from two sources. First, 
the richness and precision of the available experimental data on many observables which 
to compare with, is certain to lead to more stringent tests of nuclear potential models 
than neutron-deuteron (nd) scattering with its much smaller and much less precise data 
base. Secondly, the necessity to include the Coulomb interaction in a way that is both 
mathematically correct and practical has been, and still is, one of the outstanding theoretical 
tasks. 

Of the several approaches that have been proposed to take into account Coulomb inter- 
actions in charged-composite particle reactions, only two have reached the status to permit 
concrete numerical calculations. The most obvious one, namely to work with the Schrodinger 
(see Refs. JT|-^| and references therein) or equivalently differential Faddeev-Merkuriev equa- 
tions 0-0] in coordinate space, requires knowledge of the complete boundary conditions, in 
order to guarantee uniqueness of the solution. Below the three-body threshold when only 
two-cluster channels are open, this presents no difficulty However, above that threshold the 
complete boundary condition to be imposed in the region where all three particles even- 
tually become asymptotically free PHTT|1 has, to our knowledge, not yet been implemented 
satisfactorily in any solution scheme. 

Based on momentum-space three-body Faddeev \Y1\ or Alt-Grafiberger-Sandhas (AGS) 
|T3|j equations, mathematically well defined integral equations for charged-composite particle 
reactions have been derived for two cases. If one of the three particles is uncharged and the 
others have charges of equal sign, the AGS equations for the three-body transition operators 
have been proven by Alt, Sandhas, and Ziegelmann |14| to possess compact kernels for all 
energies in a special class of functions (see also ||). That is, these equations are amenable to 
stardard solution methods, thereby yielding the physical amplitudes for all reactions which 
are possible in such a system, at all energies. If, however, all three particles are charged (with 
charges of equal sign) only somewhat limited information is available as yet. Indeed, if the 
three-body energy is negative, i.e. below the three-body threshold, the Faddeev equations 
have been reformulated in such a way that the kernels of the new equations possess compact 
kernels | ]T5| , [T5| ] . For positive energies compactness could be proven only for the kernels of 



certain integral equations for effective-two-body transition amplitudes, i.e., for amplitudes 
which describe all possible binary (i.e., (in-)elastic and rearrangement collisions, or quite 
generally so-called 2^2) reactions flT7| , [T8|j (this proof holds a fortiori if only two particles 
are charged and one is neutral). The formulation of analogously modified integral equations 
with compact kernels for breakup (2 — > 3) or even 3^3 processes is still lacking. 

A practical solution method for the aforementioned equations for effective- two-body tran- 
sition amplitudes has been developed in |T^,|T^|,|2D| (see also the review PHI ). Starting from 
screened, and thus short-ranged, Coulomb potentials, the usual short-range equations are 
solved by standard methods. The physical amplitudes corresponding to unscreened Coulomb 
potentials are then recovered by numerically performing a limiting procedure in suitably 
renormalised quantities. Details can be found in |2l|j2^]. This approach has been applied 
with great success to the calculation of differential cross sections for pd elastic scattering 



22| , |23|| (see also [21| and references therein) and to five-fold differential cross sections for 
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pd breakup [p4H26| in various kinematic configurations, with due account of the Coulomb 
interaction (but employing simple models of the nuclear interaction only). 

Recently we have communicated the first successful calculation of proton-deuteron scat- 
tering observables for the Paris potential using this same approach j27|]28|1 . Here, results for 
many more energies will be presented and compared with experimental data. 

The plan of the paper is as follows. In Sect. [H] we give a brief recapitulation of the most 
important aspects of the screening and renormalisation method. Section |inj then contains 
results for differential cross sections and various polarisation observables. The final section 
contains our conclusions. In the Appendix the explicit expressions for the partial-wave 
decomposed effective potentials are collected. 

As usual we choose units such that h = c = 1. 



II. FORMALISM 

For the convenience of the reader we briefly recapitulate the basic equations |T1 



1. Notation 

Consider three distinguishable particles with masses m u , v = 1,2,3. Moreover, two of 
them, say particles 1 and 2, are supposed to be charged, with charges e\ and e2, satisfying 
e\C2 > 0. We use the standard notation: on a one-body quantity an index a characterizes the 
particle a, on a two-body quantity the pair of particles {(3 + 7), with /3, 7 ^ a, and finally 
on a three-body quantity the two-fragment partition a + (fiy) describing free particles a 
and the bound state (Pj). Throughout we work in the total center-of-mass system. Jacobi 
coordinates are introduced as follows: k a is the relative momentum between particles (3 
and 7, and \i a = mpm 1 / {m^ + m 7 ) their reduced mass; q Q denotes the relative momentum 
between particle a and the center of mass of the pair (j3j), the corresponding reduced mass 
being defined as M a = m a {mp + m 7 )/(m a + rap + m 7 ). 

The Hamiltonian of the three-body system is 

3 

H = H + V = H + J2 V »> C 1 ) 

with 

H = K 2 j2fi a + Q 2 j2M a (2) 

being the free three-body Hamiltonian. K a and Q Q are the momentum operators with 
eigenvalues k a and q Q , respectively. 

The two-body interaction in subsystem a has the general form 

= V Q S + 5 a3 V*, (3) 

where V^f is the 'short-range' (i.e., nuclear) part and 

V«(r) = e"*/* (4) 
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the Coulomb potential which for practical reasons we assume to be exponentially screened, 
with screening radius R. The corresponding T-operator is given as solution of the usual 
Lippmann-Schwinger equation with the full interaction V c 



(R) 
a , 



Ti R \z a ) = + V^G {z a )fi R \z a ). (5) 

For clarity, energy-dependent two-body operators, when read in the two-particle space, are 
characterized by a hat. Moreover, z a denotes the energy in subsystem a. The fact that the 
potential is a sum of two terms, cf. Eq. @, carries over to the transition operator (|5]): 

fi R \z a ) = 5 a3 f R (z a ) + f* R {z a ). (6) 

Here, T R (z a ) the pure screened Coulomb transition operator for the pair of protons 1 and 
2. The so-called Coulomb- modified short-range transition operator T^ R (z a ) is given as 

f™(z a ) = [1 + S a3 f R (z a )G (z a )}tl R (z a )[l + 5 a3 G (z 3 )f R (z 3 )\, (7) 
il R {z a ) = V! + V!Gi R \z a )Vi. (8) 

For the particle pair (/3 + 7), the free two-body resolvent is denoted by G (z a ) = (z a — 
K 2 /2fi a )~ 1 and the full resolvent by 

& R \z a ) = (z a - K 2 j2^ a - V^Y 1 = (z a - K 2 j2n a - V! - 8 a3 V R y\ (9) 

We point out that if there exists a bound state of energy B a < in subsystem a, the 
corresponding T-operator T a R \z a ), and hence also T^ R (z a ) } must have a pole of the form 

fi R \z a ) Za ~ Ba f s a R { Za ) Za ~ Ba Y^MMY^ : (10) 

Za B a 

where \ip a ) is the appropriate bound state wave function. Generalisation to several bound 
states is obvious. 



2. Equations for the three-body arrangement operators 

( Ft) 

The AGS three-body transition operator U^g (z) which leads from a partition a + (/3, 7) 
of the three particles to a partition (3 + (0,7) is defined as solution of 

3 

Ulf(z) =5 aP G Q \z)+Y.^Tl R \z)G,{z)ul R \z) (11a) 

3 

= 5 aP G^(z) + J2U a R) (z)G (z)Tl ) R \z)6 uP . (lib) 

Here, 5 a p = 1 — 5 a p is the anti-Kronecker symbol and Go(z) — (z — Ho)' 1 the resolvent of 
the three-free particle Hamiltonian H . 

The splitting ([]) of the subsystem amplitudes induces a corresponding splitting of the 
three-body operators U^g(z). Define new operators U R f3 (z) as solutions of the same AGS 
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equations (|TTD but with only the Coulomb part of the subsystem amplitudes in the kernels, 



U R p {z) = 5 aP G,\z) + J2s a uT u R (z)Go(z)U^(z) 

v=\ 

3 

u=l 

Then U^\z) and U R p{z) are related via 

U$\z) = U R ,(z) + j^Ut{z)G Q {z)T^z)G,{z)ul R \z) 

v=\ 

3 

= U%(z) +J2u^(z)Go(z)T„ SR (z)G (z)U*(z). 



(12a) 
(12b) 



(13a) 
(13b) 



An important practical simplification arises if only two particles are charged as it happens 
in the present case. For, Eqs. (^) can be solved explicitly to yield 



U R (z) = SapGz\z) + 5 3 J^T R (z). 
As a consequence, Eqs. (|T3"D with (JUty are exact. 



(14) 



3. Physical transition amplitudes 

Let the initial channel state |'0Q)|qa) be given as the product of the bound state wave 
function \ip a ) (belonging to the binding energy B a ) of the pair (/?, 7), and the plane wave 
|q a ) describing the free motion of particle a relative to the center of mass of this pair. 
Analogously for the outgoing channel state. Then the plane-wave matrix element 

T Q ( J(q a , q^; E + iO) = (q a \^ a \U {R) (E + iO)\M^ ) (15) 

is the physical transition amplitude for screened Coulomb potentials, provided the incoming 
and the outgoing energy are related to the energy parameter E via the energy-shell relation 

E = E a = q Q 2 /2M a + B a = Efs = o[^/2M p + Bp . (16) 

In order to extract the desired amplitude pertaining to unscreened Coulomb potentials 
the on-shell amplitude T^(q a , q^; E + iO) has to be multiplied by appropriate renormali- 

sation factors Z a ^ 2 (q a ) and Z '(q'p) which are uniquely determined by the special choice 
of screening function, and the limit R — > 00 has to be performed: 

r o ^(q a) qJ;J5 + iO) = jjm (fe^gta, <fr E + iO)Z$ 2 (17) 

Details of this procedure are described below. 
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4- Special case: separable nuclear interactions 



As mentioned in the Introduction, in principle the coupled equations (|TT|), with R set 
equal to infinity, for the three-body operators U a p(z) could be solved as they stand. But 
the presence of the highly singular Coulomb T-matrix in the kernel is certain to present 
formidable numerical difficulties. We, therefore, have adopted another solution strategy. 
Namely, we use separable approximations of the original (local or nonlocal) nucleon-nucleon 
potentials. That is, we assume to be represented as a sum of separable terms, 

N a 

Va = 2^ Warn) ^a,mn(^Pan\- (18) 
m,n=l 

Here, the index m of the (nuclear) form factor \tp am ) not only characterizes the complete 
set of quantum numbers which uniquely characterizes a given state of the particle pair a 
but also enumerates the number of terms per fixed set of quantum numbers, i.e. the rank. 
It will be specified later. Note that this assumption does not represent a loss of generality 
as any given short-range potential can be approximated in a form like ([18]), to any desired 
degree of accuracy. 

Let us introduce the Coulomb- modified form factor 

\9am) = [1 + 5 a3 fX(z 3 )G (z 3 )]\tp am ) (19) 

which differs from the nuclear form factor only for the pp subsystem characterized by a = 3. 
Although \g3m) will, therefore, depend on the screening radius R this dependence will, 
however, not be indicated explicitly. Then 

T* R {z a )= \g am )Al R l n (z a )(g am \, (20) 

m,n=l 

with the elements of the matrix A^\z a ) being solutions of 

+ Krn^K„\G (Za)\9av)^L(Za). (21) 

In channels where a bound state exists (in our case the deuteron for which a 7^ 3) it 
must be ascertained that the corresponding subsystem T-matrix Ti R \ or equivalently the 
appropriate element of the matrix Ai R ', shows the pole behavior (|T0|). This is guaranteed 
if the form factor which represents the bound state and is called, say, \<p a i) = since 
«^3, satisfies 

G (B a )\g al ) = G (B a )\<p al ) = a ^ 3. (22) 
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5. Off-shell equations for transition amplitudes 



Consider the quantities 

^S >/3n (z) : = (g am \G (z)uif(z)G (z)\g Pn ). (23) 

They are effective-two-body transition operators which act only in the space spanned by the 
momentum eigenstates |q a ) between the two fragments and describe all binary (so-called 
2^2) collisions. The matrix elements (q^l^m^nC^ + zO)|q^) in which the form factors in 
the initial and final state correspond to bound states and hence satisfy a condition of the 
type (p2"D , coincide on the energy shell (|T6|) with the physical amplitudes flOp. 

It is now an easy task to derive equations for X*mf3n( z )- Sandwiching Eqs. (|i~l"D between 
(g a m\ and \gp n ) yields the coupled, multichannel, Lippmann-Schwinger-type equations 

3 N-y 

tZniz) = v2^w + E E vS^w^w^lw ( 24a ) 

7,5=1 i,j=l 
3 iV 7 

= vS^w + E E ^w^wv&w- ( 24b ) 

7,5=1 i,i=l 

The effective arrangement potentials V^g n (z) are defined as 

= 8 a p(g am I G (2) I #/3n) + SapSaaifam I ^(^^(^GoOO | 

+S a pS a3 6p3((p am I G (^)7 1 3 iJ ( 2; ) G 'o( (25) 
-EVSI, (26) 

i=0 

They are depicted in Fig. |l|. Note that on account of the definition (|T9|) the first term in Eq. 
(plf), and hence also in Fig. [I], comprises actually three different contributions, enumerated 
by i = 0, 1, 2, depending on whether both form factors are purely nuclear or either one of 
them is Coulomb-modified. 

The plane wave matrix elements of the effective propagators G^anpS 2 ) are gi ven as 



6. Angular momentum decomposition 

For various reasons we found it more convenient not to work with the isospin formalism. 
Hence we use the following angular momentum coupling for a given channel a: s@ + s 7 = S Q , 
L a + S a = 3 a , s a + 3 a = Sq,, l a + S Q = J. Here, s„ denotes the spin of particle u, L a the 
relative orbital angular momentum, S a the total spin, and 3 a the total angular momentum 



7 



of particles (3 and 7; moreover, \ a denotes the relative orbital angular momentum of particle 
a and the pair (J3j), and finally J the total angular momentum of the three-body system. 

In order to simplify the notation, in the following explicit subsystem indices on partial- 
wave projected genuine two-body quantities and channel indices on effective two-body quan- 
tities are omitted. 

Let the partial wave expansion of the short range interaction ( |18|) between particles (3 
and 7 be given as 

Kf(p,p')=47r Yl y^iP^^^^iPt, (28) 

Jot,Mj a ,S a L a ,L' a = \J a — Set I 

with 



N 



v L s : J L t( P ,p')= £ d±mi±wy4fcW 



(29) 



u a ,u'=l 



Here, = Em 1o m So ^^.^(pjxwf... N ^ a * ^ rank of the sep- 

arable expansion in the two-body channel a with fixed J a and S a , fj^J^ip) are the cor- 
responding form factors which are chosen real, and \ J T " T ,„, = X T ",? T „ are the (real) 
potential strengths. The Coulomb-modified form factors are 



,J a M, 



9lZiZ(p) 



<PL a Sl(p) +Sc3 



2tt 2 



(30) 



where T 3 R Li is the screened partial wave Coulomb T-matrix. With these definitions the plane 
wave matrix elements of the Coulomb-modified short-range T-operator take the form 



p|Tf H ( Z J|p') = 4 lr 



E 

J a ,Mj a ,S a ,L 

a 1 a 



ytl Ja (p)[y J L is a Ja (p')] 



£ gt± (p)^ R lt< Msi& V) ■ 



(31) 



The matrix elements Ai „" f, a , (z a ) are obtained as solutions of 



•> ct S a. 



After partial-wave decomposition the effective potentials with given J and parity ir for 
a transition a{f3 r y) —>■ f3( r ya) can likewise be written as a sum of five terms: 



(33) 



i=0 



where 



S 



V£IZ%^ fa Z ) = E A $fa> ^) R SX,n ,Ma, Z), (34) 

K 

li a — 1, • • • N JaSa , up — 1, • • • N j i sS i 3 . The multi-indices m Q and are defined as m a =(l a , s a , 
Eq,, J a , L a , S'q, S/3, s 7 ), and n / 3=(/ ( g, s^, E^, J^, Lg, 5^, s 7 , Sq,). Finally, for the contributions 
with i = 0, 1, and 2, k = C is a single index while for i = 3 and 4 it is a multi-index 
k = (CiC 2 f). The explicit expressions for the functions At\q a -, q'p) and Rm^ a ,n p up(q a , Qp', z) 
can be found in the Appendix. 

Thus, we have to solve the following coupled set of integral equations (E + := E + iO) 

poo (iq"q" 2 ^jvr 

-V f f T t' ^ 

x (JS + - < 2 /2M 7 )^^;^ <£; £ + ) . (35) 

The meaning of the various indices has been described above. 

In general, the physical T-matrix element T z v . s . / ou tS otl t^ a ' %'■> E +) w hich describes the 
transition from channel a, where the particle pair (/?, 7) is a deuteron state with Ji n = 1 
which is not explicitly indicated, the relative orbital momentum of between particle a and 
deuteron is k n and the total channel spin is E in , to channel j3 where the particle pair (0,7) 
is in a deuteron state (J ou t = 1) and the channel orbital angular momentum and spin are 
l out and E out , respectively, can be calculated from the solutions of Eq. (^) as 

Switching off the short-range interactions will reduce the coupled equations ( |3"5j ) to those 
pertaining to the partial-wave decomposed screened center-of-mass Coulomb scattering am- 
plitude 

T KL,z otltSout (^, q'p\ E+) = MiinU^^a - ^)j^q~ a (e 2 '^ (9Q) - l) , (37) 

where crf-(q) are the screened Coulomb phase shifts. From these two amplitudes the 
Coulomb-modified short-range T-matrix follows directly as 

T LsL,/ otl iS otl i^ Q ' Qp'' E +) = T L l s„ 1 ,/ OU iS OIl i(9 , a' q'p] E + ) - T ; ^ Sm)ZotitSoiit (g a , q'p] E + ). (38) 

The unscreening procedure is now performed in the amplitude (|38|) . We multiply 

T hSZiou^ouM^ ( l^ E +) h y the ^normalisation factors Z~^ 2 (q a ) and Z^ 2 (q' p ) and re- 
peat the calculation with increasing value of the screening radius R until the result has 
become independent of it. In this way we end up with the unscreened Coulomb-modified 
short-range amplitude 

T kSLiout^ouM^ fa E+) = Jim ^ai /2 (^)T{f i f^ /outSoiit (g a , q^ E + )Z^{q^. (39) 

Finally, summing up the partial wave series and adding to the result the analytically known 
(unscreened) center-of-mass Coulomb scattering amplitude yields the final reaction ampli- 
tudes from which the various observables can be calculated. 

We emphasize once more that solution of Eq. (^) yields, after execution of the unscreen- 
ing procedure as described above, charged-composite particle transition amplitudes which 
are exact for a nuclear potential of the form ([18]). 
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III. RESULTS 



As mentioned in the Introduction, the first theoretically satisfactory calculations of pd 
elastic scattering above the breakup threshold p3|j2^] and of pd breakup [f24]-|2"6|| employed 
rather simple ansaetze for the nuclear interaction. In spite of this limitation at least semi- 
quantitative agreement with experimental differential cross sections for elastic scattering and 
for five-fold differential cross sections for deuteron breakup in various kinematic situations 
could be achieved. However, for a more detailed comparison with experimental data, in 
particular for polarisation observables, more sophisticated nuclear potential models must be 
used. For this reason we have performed calculations with the realistic Paris potential. First 
results have been published recently |27| , |28| . Here, we present some extended calculations 
of differential cross sections and various polarisation observables for elastic pd scattering. 

We used the Paris potential in the PEST1 form ||29|| . S and P waves were included in 
the pp and the np spin singlet channels, and the coupled S-D waves in the np spin triplet 
channel. This leads to maximally 29 coupled integral equations to be solved. The number 
of total angular momenta in the pd system was chosen so high that stable results for all 
observables were obtained. It was found that J = 17/2 suffices for the lower energies, and 
J = 19/2 for the two highest energies, for the level of accuracy aspired to. 

In order to perform the unscreening of the resulting amplitudes numerically, repeated 
solution of the integral equations (p5|) with and without the nuclear interaction is required, 
with the screening radius R increased until the r. h. side of Eq. ( P5| ) becomes independent 
of it. We found that R = 625 fm was enough for all purposes. 

When calculating the effective potentials (|4]) and effective free Green functions ([32]) 
we have made only one approximation. Namely, as indicated in the explicit expressions 
given in the Appendix, we have used the Born approximation for the pp Coulomb T-matrix. 
As has been shown in when the range of the form factors is of typical nuclear size, 
this approximation is accurate to a few percent for all energies and scattering angles (in 
contrast to the atomic case where it typically fails by several orders of magnitude). And 
since the Coulomb interaction modifies the purely nuclear pd phase parameters by at most 
10 percent, the error introduced by this approximation is therefore estimated to be well 
below the 1 percent level. 

In Fig. [2] we present differential cross sections for proton laboratory energies from 2.5 to 
22.7 MeV. For comparison the corresponding results for nd scattering are included. Inspec- 



tion reveals that very good agreement with the experimental data of Sagara et al. |3l[] and 



of Sperison et al. |32] is achieved, except at the lower energies where our calculations slightly 
overestimate the data. For the vector analyzing powers depicted in Fig. |3] the reproduction 
of the data is much less satisfactory. In particular, at the lower energies the maximum of 
the vector analyzing power is strongly underestimated. This is the so-called the A^-puzzle 
which has been with us for a long time in neutron-deuteron scattering and is also present 
in the pd reaction as already noted in [0] . In spite of a variety of speculations regarding its 
origin and remedy (see, e.g., p?5]-|3"7[ and references therein), at present no satisfactory solu- 
tion to this problem is available. But it appears that the failure of the theory to reproduce 
the experimental maximum disappears at higher energies, at the expense of an increasing 
discrepancy in the minimum around 100 degrees. A similar situation occurs for iTn as can 
be inferred from Fig. [|. Experimental tensor polarisations T 20 ,T 2 i, and T 22 where available 
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are reasonably well reproduced by our calculations as can be seen from Figs. | - 0. For all 
observables presented the modifications due to the Coulomb interaction are rather strong 
at the lower energies but eventually become small although not negligible at the highest 
energy. We mention that our results are rather close to those of Ref. @]. 



In Ref. |31] it was pointed out that apparently the magnitude of the experimental dif- 
ferential cross section minimum differed appreciably from theoretical results. In fact, the 
relative difference A min := (o"™ e ° r — cr™™)/^™ was found to be rather large and positive at 
low energies, to change sign around 5 MeV and to become negative large at higher energies, 
reaching - 25% at 18 MeV. Explanation of such a strong, and strongly energy-dependent, 
effect seemed to be very difficult. It was therefore suggested that there exists another real 
discrepancy, later termed 'Sagara discrepancy', between experiment and theory, besides the 
Ay-puzzle. This inference was, however, not very compelling as the theoretical calculations 
used for comparison had actually been performed for nd scattering, with only a very rough 
account of Coulomb effects. The existence of such an energy-dependent discrepancy was 
later corroborated by calculations employing several 'realistic' nuclear potentials but again 
relying on the same approximation for including Coulomb effects |[33|| . Thus, in both calcula- 
tions it was ignored that this Coulomb correction method had already been demonstrated in 



23J to be generally unsatisfactory and, in addition, to lead to a rather strong energy depen- 



dence of the failure, in particular in the cross section minimum. A further attempt [[34] to 
explain this effect by A-isobar induced three-body forces ignored Coulomb effects altogether 
(but gave sizeable corrections particularly for higher energies not considered here). 

In Fig. |8| we present the relative difference A m i n for the cross section minimum between 
our theoretical results and the data of Ref. [ 31] . Inspection reveals that even with a correct 



description of the Coulomb repulsion between the protons the 'Sagara discrepancy' survives, 
albeit with greatly reduced overall magnitude as compared to the calculations with improper 
account of the Coulomb interaction. In addition, the percentage excess in A m i n has become 
only rather weakly dependent on energy in the range considered (from 9.3% at 5 MeV 
to 6.3% at 18 MeV), in contrast to the original estimates |j3l . It is interesting to note 



that for all energies considered, our calculations yield a larger cross section minimum than 
experiment, i.e., A min > 0. However, before drawing any conclusions about the origin of 
this overestimation it should be kept in mind that in particular the cross section minimum 
is very sensitive to the finer details of the nuclear force model. Hence, it could well be that 
in calculations using higher-rank, and thus better, approximations of the Paris (or a more 
modern) potential, the remaining difference even disappears. 



IV. SUMMARY 

In this work we have presented the results of our calculation of differential cross sec- 
tions and polarisation observables for proton-deuteron elastic scattering in the energy range 
from 2.5 to 22.7 MeV, with due allowance of the Coulomb repulsion between the protons. 
The mathematical framework for rigorously taking into account the long-ranged Coulomb 
interaction in momentum space integral equations is provided both by the screening and 
renormalization approach [TJ]] and by the investigation of the analyticity properties of the 



kernels of the pertinent equations |17|JT8[] . The former even provides for a practical solution 
scheme which has been adopted in the present work. 
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As input for the nuclear interaction we employed the PEST1 version of the Paris poten- 
tial which is well known to represent an excellent (separable) approximation to the original 
potential. The calculated differential cross sections led to a very satisfactory reproduction 
of the experimental data. This fact also gives rise to a decisive reduction in absolute magni- 
tude and in its energy dependence, of the so-called 'Sagara discrepancy', originally described 
in Ref. ||31|| . It even suggests that the latter might cease to exist when more sophisticated 
nuclear potential models will be used, provided due account is made of the Coulomb inter- 
action. For the various vector and tensor polarisation observables the agreement is not as 
good, as was to be expected from the fact that a similar lack of agreement is known to occur 
in neutron-deuteron scattering. 

In order to shed some light on the origin of the remaining discrepancies, calculations 
with improved nuclear input are called for. Such are under way. 



APPENDIX: 



In the Appendix we present the explicit expressions for the various contributions to 
the angular momentum projected effective potential. The following notations are used: 
[/] = \j2l + 1, [I 2 ] = 21 + 1, X a p = m a /(m a + mp) = 1 - Xp a , a ^ (3. Moreover, e a p = -ep a 
is the antisymmetric symbol with e a p = +1 if (a, 0) is a cyclic ordering of the indices (1,2,3). 
Unit vectors are denoted by a hat, k = k//c. 



A 



(0,1,2) 



^Sct—sp—s-f 5/3— s 



z af3 



[l a L a S a Yj a J a lpLpSpYjpJpJ 2 £ 2] 

47T 



X 



IMA. - Aa)A^(^ - A ))qT^qT " " l "A:: A 



A a +A /3 ^/(-L Q +L / 3)-(A a +A (3 ) x \ a i^-Aj 



X, 



(2^ + 1)1(2^ + 1)! 



(2A a + l)!(2(L a - Aq,) + l)!(2A /3 + l)!(2(L /3 - A^) + 1)! 

r2 f 21 



AhM 2 f 



X 







M 2 M 1 f 

la 1/3 £ 




A a Ap Mt 




Mi l Q C 




A a A a 



Lp Lp — Ap Ap 
f M 2 M 1 



M 2 lp C 


S a s a Sp Lp 
J a S-y Jp f 
L a S a sp £/? 



1 f + 

Rt!S C n u ((la,q'p; Z) = 5 a p- J ^ dxP C [ 



X 



X 



k^kp L "git(k a )g^(k'p) 
z - q 2 a /2M a - k 2 j2^ a 



(Al) 



12 



with k a = e a p(X 01 q a + q^), k' = e /37 (A Cf7 q^ + q a ), x = q a • q^. 



/|( 3 ) ( n n > \ — X X L a +L fj( ^j_ Sa+Sa+ l a+ l g+ C 2 +f [IgEaJJp'E'pJpJ j£?A 



X- 



(2^)! f Z Q S Q 



(2(L /3 - /))! i ^ l p f j { Jfs Hp f j { L a L p S a 



x 



^ y/(2(f-w))\(2w)\ 



W 1/3 £ 2 \ ( L a d Lp - f\ j f - w l a £ 2 \ ( L f - w C 2 
J \ J \ h w f j \ 

v V" r L P ML f3 r L P M Lp r L a M La fA , 

X ^fM Lp -M La L a M L J^Lp-fM La fM L0 -M Lc yLp-fM La CxO-i ^) 

M La M Lg 

with 

r L af iLpuf {Ha, q/J, - 



dkk 2+L ?-fg J L f^(k) /•+idx 1 P £l (x 1 )^(|A a + k|)|A Q + k'- i 



o 2 - g2/2M Q - fc 2 /2// Q J_ x z - qg/2M p - (A a + k) 2 /2^ 



(A4) 



Here, A^ = q a - q^, A a = -A 7/3 A^, a; 2 = q a • q'p, £i = q ■ A a . 

In order to increase the numerical accuracy it is advantageous to extract, in the deuteron 
channel, from expression (|A3| ) the pure screened center-of-mass Coulomb potential. For this 
purpose we only have to take into account that for n = m and q a = q^, i.e., A a = A' a = 0, 

the only non-vanishing element of expression (|A2|) is Aj^ C2Q (q a7 q'g) = 5 £ l0 <5 c 2 i a I ■ Hence, 



in the deuteron channel we introduce for q^ = q~ the short-hand notation 

Qa! z) = 16tt 3 F JaSa L af ,u(q a ; z), (A5) 

with 

i °2 k 2 dkg J T aSa (k)g J T aSa (k) 

I (z -A:l k t%J)* deuteron channel, 



Fj a Sa!* ia ,(qa\z) = { {2n)3J o (*-<&/™*-kW — , (A6) 

elsewhere. 

The requirement that the deuteron wave function be normalized to unity implies 
J2L a nu Fj a s a L atJ ,u(qa', E a ) = 1. Consequently, expression ([AID for can be rewritten 

as 
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+ 



57T 



J^ 1 dx 2 p C2 (x 2 )vi R \A' a ) (FtS?(q„q^) - ^ 3 F JaSaLa ^z] 



(A7) 



Here, is the partial-wave projection of the screened Coulomb potential. Thus, Eq. (|A7|) 
is in a numerically convenient form: the screened center-of-mass Coulomb potential which in 
the limit R — > oo gives rise to the troublesome 'Coulomb singularity' is explicitly extracted 
in the simple first term, while in the second term which in general has to be computed 
numerically the latter does not longer occur. As an additional bonus the decomposition 
flff7|) provides us with the possibility to use the same set of coupled equations fl55|) to find 
screened center-of-mass pure Coulomb scattering amplitude. One simply has to switch off 
all nuclear interactions, keep as effective potential only the first term of ( |A7| ), and replace 
the effective free Green function A^jl, S J,(E + — q" 2 /2M~ f ) by the free 'proton-point deuteron 
propagator' 1/(E + — B d — g" 2 /2M 7 ), where B d is the deuteron binding energy. 
Finally, 



.4 



(i) 



La+SaSp-S-t Lg+Sg-Sa-S^ /_ 1 n. j + s l3 +s l3 +2s~ l +l a +L l3 +l l3 +C2+f 
■ya t/5„, \ J- I 



[laYjaJaSalpTipJpSpLlfCl] \J {2Lp)\ 



X ■ 



47T 

^a S a S/3 Lp 



(2(^-/))! 



J, 



a ^7 J/3 f 
L a S a S/3 Tj/3 




X 



^ y/(2(f-w))\(2w)\ 



1 w 


b 


^ ( 


V o 








L a £l L/3 — f 





/ - W l a C 2 

h W f 



l a f ~W C 2 





L P M Lg 



C 



LaM L , 



c 



L a M La 



fM L -M La L a M La ^ Lg-fM La fM L -M La ^ Lg-fM^dO 



(A8) 



d(4)£i£ 2 / („ „' 



+i 



dx 2 Pc 2 (x 2 ) 



00 rlh h 2+Ll3 ~f n JaSa (h 



X 



z - q 2 a /2M a - k 2 a /2^ c 
+i d Xl P Cl (xM R Kk' a , k a , x 1? • q)g J L ff(\k a + q|)|k a + q|~^ 



-i 



(A9) 



z - q^/2M - (k a + q) 2 /2/i/3 

where q = (-A 7/3 q a + A 7a q^), k' a = (q^ + A^q^), x 2 = q Q ■ q' p , x 1 = k a ■ q and 
v( R \k',k,x,y) = l/v 7 ^' 2 + k 2 + 2k'kxy + 1/i? 2 ) 2 - Ak' 2 k 2 {l - x 2 ){l - y 2 ). 
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FIGURES 




FIG. 1. Graphical representation of the effective potential (|26|). The first diagram represents 
the three terms Vma^navai depending on whether both form factors are purely nuclear (open 
semicircles) or either one of them is Coulomb-modified. 
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FIG. 2. Proton-deuteron (solid lines) differential cross section for several projectile energies. 
Experimental data are from Refs. and |32| . For comparison also the results for neutron-deuteron 
scattering (dashed lines) are given. 



18 



0.08 : 


2.5 MeV 


5 MeV 


6 MeV / - 


0.06 : 
















0.04 : 




J\ - 


,i j V 


0.02 : 


'/ \* 1 

/ V 

■"*/ \- 




\ " 




8 MeV ; 


10 MeV : 


12 MeV ( m ; 


0.15 ; 




ft 


7 V ' 


< 0.10 - 




V W : 








/ \ ' 


0.05 - 




7 \ ^ 






14 MeV ; 


18 MeV ; , ; 


22.7 MeV , : 


0.2 : 


« 

/ V 




h 


<T 0.1 : 








0.0 ; 






, lit 




60 120 


60 120 


60 120 




®c.m. 


®c.m. 


®c.m. 



FIG. 3. Nucleon vector analyzing power A y . Notation as in Fig. 
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FIG. 8. 'Sagara discrepancy' as function of the proton laboratory energy. Open squares: Ref. 
1[, black squares: present calculation. 
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